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Abstract 

Practical Experiences witli Films, Cameras and 
Navigation Systems 

R. W. Graham & R. E. Read 
^Aerial Imaging Systems Ltd. 
Abstract 


The first author commenced his experiences with aerial 
photography in 1945, when he joined the Royal Air Force (RAF). 
The aerial cameras used at that time were the Williamson F24 
and F52, with Ilford HP3 and FP3 panchromatic films. 
Subsequently the Williamson F95 camera was tested against the 
USAF Sonne 11 camera. During his final years in the RAF, a 
variety of film types was employed: colour, monochrome infra- 
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red and false colour infra-red. Additional experience was gained 
when the author was employed by the International Training 
Centre (ITC), The Netherlands and he saw the further 
development of the air camera from the Wild RC8 and Zeiss 
(Oberkochen) RMK 15/23 to the Wild RC20, Zeiss (Jena) LMK 
and Zeiss (Oberkochen) RMK TOP. He predicts that the future 
lies with digital cameras. 

The second author started working with aerial photography in 
the 1950s. Details are given of cameras, visual navigation sights 
and survey aircraft dating from that period. Overseas operations, 
carried out mainly by Hunting Surveys and Fairey Surveys, 
posed particular problems because the majority of the work 
there was undertaken without the aid of existing maps. The 
introduction of external navigation systems is described, from 
the Decca ship navigation system, through the GNS-200 
VLF/Omega low frequency radio system to INS, the Inertlal 
Navigation System. The author later joined the ITC, which 
enabled him to participate in trials of the Computer-controlled 
Photo Navigation System. Eventually practical experience was 
gained with differential GPS, particularly in countries in 
southeast Asia. Aerial photography currently appears to benefit 
less from financial investment than the photographic, aviation 
and survey industries. However, the satisfaction of making a 
significant contribution to major overseas development projects 
outweighs the frustration of lack of acceptance of plans for the 
future of the air survey industry. 
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What is a line scan camera? 


A linejcanxamejajs an image capturing device having a CCD sensor which is 
formed by a single line of photosensitive elements (pixels). Therefore, unlike area 
sensors which generate frames, in this case the image acquisition is made line by line. 

One single scanning line of a line scan device can be considered as a one-dimensional 
mapping of the brightness related to every single point of an observed line. A linear 
scanning generates a line, showing on the Y axis the brightness of each point given in 
grey levels (from 0 to 255 levels), 

A sudden change of the grey level in a single point corresponds either to a point on 
the edge of an object or to any colour or aspect variation of the acquired image. 
Detection of this change allows a precision measurement, thanks also to the high 
resolution on the linear sensor which is considerably better than the resolution of an 
area imager. 


For instance, by using a backlight you may know the position and width of a strip, by 
observing the transition of every single point; or detect the impurities (holes, 
scratches, spots etc.) on any inspected surface. 

For this type of camera the use of lenses for line scan cameras are reconmiended. 
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Abstract: 

A computerized flight inspection system is disclosed. The system 
of the present invention may be utilized to generate an accurate 
reference location with respect to an airport runway for an aircraft 
having an inertial navigation system. A selected geometric pattem 
having a highly unambiguous autocorrelation function is placed on 
at least one end of the runway, ^video line scanning camera ^ 
mounted to the aircraft is then utilized Wscain the geometric 
pattem in a line generally perpendicular^^ the line of flight. The 
output of the video li ne sca n n ingxamera is then correlated with a 
stored reference pattern to generate a signaFindicative of the 
detection of the geometric pattem on the runway. A laser altimeter 
is mounted to the aircraft and utilized to generate an accurate 
signal indicative of the aircraft altitude with respect to the runway 
pattem. The outputs of the correlation circuit and the laser 
altimeter are then utilized to correct data from the inertial 
navigation system. 
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This is a continuation-in-part of co-pending application Ser. No. 
065,075, filed on 6-17-87, abandoned, which is a continuation-in- 
part of application Ser. No. 946,124, filed 12-24-86, abandoned, 
which is a continuation-in-part of the present application Ser. No. 
699,333 filed 2-7-85, abandoned. 

What is claimed is: 

1 . A computerized flight inspection system for providing an 
accurate reference location with respect to a runway for an aircraft 
having an inertial navigation system, said flight inspection system 
comprising: 


selected indicia disposed on at least one end of an airport runway; 
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video scanning means disposed within said aircraft for scanning 
said selected indicia in a line generally perpendicular to said 
aircraft's line of flight and for generating an output indicative of 
the scanned video pattern; 

memory means for storing a reference video pattern; 

correlator means coupled to said video scanning means and said 
memory means for correlating said scanned video pattem and said 
reference video pattem; 

altimeter means for measuring the altitude of said aircraft and for 
generating a signal indicative thereof; and 

processor means coupled to said correlator means and said 
altimeter means for generating correction data for said inertial 
navigation system in response to said correlation and said altitude 
measurement. 

2. The computerized flight inspection system according to claim 1 
wherein said selected indicia comprises a selected geometric 
pattem. 

3. The computerized flight inspection system according to claim 1 
wherein said video scanning means comprises a solid state line 
scan camera. 

4. The computerized flight inspection system according to claim 3 
wherein said solid state line scan camera includes a charge couple 
device line scan image sensor. 

5. The computerized flight inspection system according to claim 4 
wherein said charge couple device line scan image sensor includes 
a single row of 1,024 image sensing elements. 

6. The computerized flight inspection system according to claim 1 
wherein said correlator means comprises a sixty-four bit 
monolithic digital correlator. 

7. The computerized flight inspection system according to claim 1 
wherein said altimeter means comprises a laser altimeter. 

8. The computerized flight inspection system according to claim 1 
wherein said processor means comprises an appropriately 
programmed microprocessor. 

9. The computerized flight inspection system according to claim 1 
wherein said memory means comprises a parallel-in-serial-out shift 
register. 
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10. The computerized flight inspection system according to claim 

I wherein said selected indicia disposed on at least one end of an 
airport runway comprises two identical indicia disposed on 
opposite ends of an airport runway. 

1 1 . A computerized flight inspection system for providing an 
accurate reference location with respect to a runway for an aircraft 
having an inertial navigation system, said flight inspection system 
comprising: 

selected indicia disposed on at least one end of an airport runway; 

video scanning means disposed within said aircraft for scanning 
said selected indicia in a Hne generally perpendicular to said 
aircraft's line of flight and for generating an output indicative of 
the scanned video pattem; 

inertial navigation means for generating a plurality of signals 
indicative of aircraft position and attitude; 

memory means for storing a reference video pattem; 

correlator means coupled to said video scarming means and said 
memory means for correlating said scaimed video pattem and said 
reference video pattem; 

altimeter means for measuring the altitude of said aircraft and for 
generating a signal indicative thereof; and 

processor means coupled to said inertial navigation system, said 
memory means, said correlator means and said altimeter means for 
modifying said reference video pattem in response to said plurality 
of signals indicative of aircraft position and attitude and said 
altitude signal; and for generating correction data for said inertial 
navigation means in response to the altitude and attitude of said 
aircraft and the output of said correlator means. 

12. The computerized flight inspection system according to claim 

I I wherein said selected indicia comprises a selected geometric 
pattem. 

13. The computerized flight inspection system according to claim 
1 1 wherein said video scan ning m eans comprises a solid state line 
scan camera. 4ps^""'^'"" 


14. The computerized flight inspection system according to claim 

13 wherein said solid state line scan camera includes a charge 
couple device line scan image sensor. 

15. The computerized flight inspection system according to claim 

14 wherein said charge couple device line scan image sensor 
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includes a single row of 1,024 image sensing elements. 

16. The computerized flight inspection system according to claim 
1 1 wherein said correlator means comprises a sixty-four bit 
monolithic digital correlator. 

17. The computerized flight inspection system according to claim 
1 1 wherein said altimeter means comprises a laser altimeter. 

18. The computerized flight inspection system according to claim 
1 1 wherein said processor means comprises an appropriately 
programmed microprocessor. 

19. The computerized flight inspection system according to claim 
1 1 wherein said memory means comprises a parallel-in-serial-out 
shift register. 

20. The computerized flight inspection system accordmg to claim 
1 1 wherein said selected indicia disposed on at least one end of an 
airport runway comprises two mirror image indicia disposed on 
opposite ends of an airport runway. 

Description: 

BACKGROUND OF THE INVENTION 

The present invention relates in general to aircraft navigation 
systems and in particular to systems for generating an accurate 
reference location with a respect to a runway. Still more 
particularly, the present invention relates to systems for generating 
position correction signals from optical observation for utilization 
in the inspection of instrument landing systems. 

Instrument landing systems are well known in the prior art. Such 
systems can be utilized to permit appropriately equipped aircraft to 
land on airport runways under highly hazardous conditions. 
However, in order to utilize such systems it is necessary to 
periodically inspect and calibrate these systems to ensure that the 
electronic indications of aircraft location are sufficiently accurate 
vnih respect to actual aircraft location. 

Highly accurate inertial navigation systems are known which 
generate relatively accurate indications of aircraft location from a 
selected point; however, these systems are not sufficiently accurate 
to permit the evaluation of Instrument Landing Systems (ILS). 
Thus, a more accurate method of generating aircraft location must 
be utilized. 

In order to generate highly accurate reference locations of 
inspection aircraft, prior art systems have utilized a variety of 
approaches. For example, one system utilizes a ground-based 
optical tracking device such as a theodolite to generate a reference 
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location which may be compared to the electronically generated 
location provided by the instrument landing system. Other systems 
utilize ground-based optical reflectors which can be utilized to 
allow an aircraft equipped with light sources and detectors to 
obtain a reference location. 

All such known systems require extensive groundbased equipment 
and/or operators to generate a reference location of the inspection 
aircraft. As such, some of these systems are subject to human error 
and some can only be operated during daylight hours in acceptable 
weather conditions. Thus, it should be apparent that a need exists 
for a flight inspection system which does not require extensive 
ground-based equipment or personnel and which generates a 
highly accurate reference location with respect to an airport 
runway. 

SUMMARY OF THE INVENTION 

It is therefore one object of the present invention to provide an 
improved flight inspection system. 

It is another object of the present invention to provide an improved 
flight inspection system which does not require extensive ground- 
based equipment. 

It is yet another object of the present invention to provide an 
improved flight inspection system which does not require groimd- 
based support personnel. 

It is still another object of the present invention to provide an 
improved flight inspection system which utilizes optical 
observations to generate position correction signals for utilization 
in the inspection of instrument landing systems. 

The foregoing objects are achieved as is now described. A selected 
geometric pattern having a highly unambiguous autocorrelation 
fimction ij^placed on at least one end of an airport runway. A video 
d-iftei^eanmn camera jnpunted to the aircraft is then utilized to 
^can the geometf ic pattern in a line generally perpendicular to the 
line of flight. The output of the video line scanning camera is then 
correlated with a stored reference pattern to generate a signal 
indicative of the detection of the geometric pattern on the runway. 
A laser altimeter is mounted to the aircraft and utilized to generate 
^Sif^^iErate signal indicative of the aircraft altitude with respect to 
f the runway pattern. . 

The outputs of the correlation circuit and the laser altimeter are 
then utilized to correct data firom the inertial navigation system. 

BRIEF DESCRIPTION OF THE DRAWINGS 
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The novel features believed characteristic of the invention are set 
forth in the appended claims. The invention itself; however, as well 
as a preferred mode of use, further objects and advantages thereof, 
will best be understood by reference to the following detailed 
description of an illustrative embodiment when read in conjunction 
with the accompanying drawings, wherein: 

FIG. 1 is a perspective view of an airport runway equipped with 
the novel flight inspection feature of the present invention; 

FIG. 2 is a side view of an aircraft equipped with the novel flight 
inspection system of the present invention; 

FIG. 3 is a plan view of a runway pattern which may be utilized 
with the novel flight inspection system of the present invention; 

FIG. 4 is a diagram of the autocorrelation function of the runway 
pattern depicted in FIG. 3; 

FIG. 5 is a block diagram of the novel flight inspection system of 
the present invention; 

FIG. 6 is a block diagram of the video position references system 
of the novel flight inspection system of the present invention; 

FIG. 7 is a more detailed block diagram of the video position 
reference system of the novel flight inspection system of the 
present invention; 

FIG. 8 is a diagrammatic end view of an aircraft and runway 
having the runway pattern depicted in FIG. 3; 

FIG. 9 is a diagrammatic side view of an aircraft and runway 
having the runway pattern depicted in FIG. 3; and 

FIG. 10 is a flow chart illustrating the steps for using the optical 
pattern recognition system position coordinates to update an 
inertial navigation system error estimate. 

DETAILED DESCRIPTION OF THE PREFERRED 
EMBODIMENT 

With reference now to the figures and in particular with reference 
to FIG. 1 there is depicted a perspective view of an airport runway 
equipped with the novel flight inspection feature of the present 
invention. As can be seen, runway 10 is depicted as having a 
geometric pattem 12 disposed on each end thereof. As will be 
explained in greater detail herein, geometric pattem 12 is selected 
to have a highly imambiguous autocorrelation function and must 
be positioned across runway 10 with an accuracy of 0.5 inch for 
Category III system inspections. Geometric pattem 12 is preferably 
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located at the extreme ends of runway 10 to reduce wear due to 
aircraft wheel marks and exhaust gases. Geometric pattern 12 must 
not be obscured by dust, sand or snow in order to operate the 
system of the present invention; however, geometric pattern 12 is 
preferably approximately six feet wide and eighteen feet long and 
thus the area to be kept clear is quite small Additionally, because 
of the natiire of the autocorrelation operation, the system is tolerant 
to a considerable degree of pattem degradation. 

Aircraft 14 is utilized to fly over geometric pattem 12 to obtain 
accurate reference position data. Preferably, the novel system of 
the present invention is initiated as aircraft 14 enters Instrument 
Landing System (ILS) envelope 16 by middle marker 18. Also 
depicted is inner marker 20, as those skilled in the art will 
appreciate. 

Referring now to FIG. 2, there is depicted a side view of an aircraft 
14 which illustrates certain features of the novel flight inspection 
system of the present invention. As can be seen, as aircraft 14 
passes over runway 10 at an altitude of approximately thirty-five to 
seventy feet, a laser altimeter 22 is utilized to measure the altitude 
of aircraft 14. Those skilled in the art will appreciate that such 
laser altimeters are readily available in the industry and will 
provide altitude accuracy within a few inches. Laser altimeter 22 is 
utilized as an input to the positiondetermining computations of the 
present invention. Also depicted in FIG. 2 is digital scanning 
camera 24 which is utilized to scan geometric pattem 12 to provide 
reference position data as will be explained herein. As an 
installation consideration, scanning camera 24 should be 
positioned as closely spaced as possible to aircraft 14's inertial 
navigation system. As will be explained herein, roll, pitch and yaw 
infomiation firom the inertial navigation system will be utilized to 
correct the output of line scanning camera 24 and such an 
installation will reduce errors due to fiiselage deflections. 

With reference now to FIG. 3, there is depicted a plan view of 
geometric pattem 12 as installed on runway 10. While a specific 
pattem is depicted, the particular pattem selected is merely one of 
many patterns which exhibit a highly unambiguous autocorrelation 
fimction. The pattem depicted is constructed utilizing bright white 
bands 26, 28, 30 and 32, which are separated by black bands 34, 36 
and 38. In the particular pattem depicted, band 26 is three units 
wide, band 28 is one unit wide and bands 30 and 32 are each two 
units wide. Similarly, black band 34 is two units wide and bands 
36 and 38 are each one unit wide. It should also be noted that the 
aforementioned bands are not bounded by straight lines, but rather 
three separate contiguous segments of straight lines which change 
their rate of convergence at intersecting lines 40 and 42 to 
approximate an exponential function. 

Referring now to FIG. 4, there is depicted a diagram of the 
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autocorrelation function of geometric pattem 12 of FIG. 3. As can 
be seen, peak 44 of the autocorrelation function is almost twice as 
great in magnitude than all remaining peaks in the autocorrelation 
function. In this maimer, a reference pattem for geometric pattem 
12 can be correlated with a scan across pattem 12 which is 
generally parallel to lines 40 and 42 and the resultant correlation 
function will indicate with a high degree of lateral resolution and 
certainty when a "match" has occurred. 

With reference now to FIG. 5, there is depicted a block diagram of 
the novel flight inspection system of the present invention. As can 
be seen, the flight inspection system is controlled by central 
processor 46. Central processor 46 is preferably a suitably 
programmed single-band computer, such as the model HK-68-3, 
manufactured by Heurikon Corporation of Madison, Wis. 

Coupled to central processor 46 are laser altimeter 22 and precision 
distance measuring equipment (DME) 48 which are utilized by 
central processor 46, together with the output of video position 
reference system 52 to generate an accurate position with respect 
to runway 10 which may be utilized to correct Inertial Navigation 
System (INS) 50 for position offset and drift. In this manner, 
aircraft position data generated by INS 50 is thus made sufficiently 
accurate to permit evaluation of Category III Instrument. 

In a preferred mode of the present invention, INS 50 is preferably a 
Delco Carousel IV Inertial Navigation System. The Delco 
Carousel IV INS provides approximately an order of magnitude of 
improvement over conventional INS systems with regard to zero 
offset of the "x" axis and the "y" axis accelerameters. Also coupled 
to central processor 46 are strip chart recorder 58 and display 60. 
Strip chart recorder 58 may be implemented utilizing any suitable 
recorder and may be utilized to provide complete recording of each 
test flight. 

The use of a wide range of avionics equipment 62 is 
acconmiodated with this preferred mode of the present invention. 
Test receivers/transmitters for the checking of ILS (Instrument 
Landing Systems), VOR (VHF Omnidirection Range), TACAN, 
MKR (Marker Beacons), and conraiunication equipment are 
connected to an Inspection Data Acquisition subsystem 56 
containing suitable translation interfaces and a slave 
microprocessor which organizes, formats and controls the data for 
subsequent use by the central processor 46. 

All avionics data received by the central processor, plus positional 
data computed thereby, are available for recording on a magnetics 
tape unit 54 for long-term storage. These data are also available for 
visual presentation on the CRT display unit 60 or strip chart 
recorder 58. 
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Referring now to FIG. 6, there is depicted a block diagram of video 
position reference system 52 of FIG. 5. As can be seen video 
signals indicative of the scanned pattern are generated by line 
scanning camera 24. These signals are then processed by video 
processor 64. Video processor 64 is utilized to produce a quantized 
serial data stream having a bit state indicative of the detection of 
"black" or "white" by digital Une scanning camera 24. 

The serial data stream from video processor 64 are then coupled to 
correlator 66 and compared to a stored reference pattern therein. In 
a preferred mode of the present invention, correlator 66 is provided 
by a monolithic 64 bit digital correlator such as the TDC 1023 J, 
manufactured by TRW, Inc. of Redondo Beach, Calif A 
correlation score is generated for each comparison and that 
correlation score is then compared to the highest correlation score 
previously obtained and stored in digital peak detector 68. If the 
current correlation score is higher than the score stored in digital 
peak detector 68, the new score is stored within digital peak 
detector 68 and strobe signals are generated which update the data 
in both cross-track register 70 and scan number register 72. 

As each strobe signal is coupled to cross-track register 70, that 
register stores the pixel count generated by timing generator 74. 
Timing generator 74 generates a pixel count beginning with zero at 
the beginning of each scan of digital line scanning camera 24. 
Thus, the value stored in cross-track register 70 will correspond to 
the lateral position of the runway reference pattern within each 
scan of digital scanning camera 24. 

Similarly, scan number register 72 acquires and holds the number 
of the scan during which the maximum correlation score occurs. 
The scan register is generated by timing generator 74 and is set to 
zero each time central processor 46 reads a frame of data, and 
advances at the scan rate. In a preferred embodiment of the present 
invention, the scan rate is preferably 200 scans per second. 
Additionally, digital peak detector 68 is also reset to zero at the 
beginning of each data frame. 

The occurrence of the runway pattern detection is selectively 
defined as the highest correlation score above some predetermined 
threshold value. For example, a correlation score of 56 or more bit 
matches out of 64. Such a threshold will be high enough to 
discriminate against random background image pattems and will 
also tolerate some degree of reference pattern degradation. An 
event detection permits the pattem lateral position (determined by 
the contents of cross-track register 70) and the time of occurrence 
(determined by the content of scan number register 72) to be 
combined by central processor 46 with aircraft altitude (as detected 
by laser altimeter 22) and aircraft attitude (as detected by Inertial 
Navigation System 50) to compute the actual position of the 
aircraft relative to runway 10 with a precision sufficient for 


http://www.freepatentsonline.com/4792904.html 


5/24/07 


Computerized flight inspection system - Patent 4792904 


Page 12 of 16 


instrument landing system inspection. 

With reference now to FIGS. 8 and 9 an aircraft and runway 
having the reference pattern template of the present invention is 
depicted to explain the computation of aircraft position relative to 
the runway. The position of the flight-inspection aircraft is defined 
with respect to a reference point on the runway where the 
centerline of the runway intersects the threshold line. The position 
of the aircraft is described in Cartesian coordinates as follows: 

x=displacement from the centerline; 

y=displacement from the threshold along the centerline; 

z=altitude above the runway surface. 

The runway surface is assumed for calculations to be essentially 
flat and level. During flight inspection path profiles the aircraft 
pitch angle is relatively small, less than 60 and the aircraft roll 
angle is also relatively small, less than lO.degree.. Therefore, 
secondorder coupling terms are negligibly small. By applying 
small-angle approximations, the altitude above the runway, with 
precision sufficient for instnmient landing system inspection, is 
given by: 

z=T cos (.alpha.. sup.2 +.beta..sup.2).sup.l/2 
where: 

r=range to the runway surface measured by the laser altimeter; 
.alpha.=aircraft roll angle; 
.beta.=aircraft pitch angle. 

By applying small-angle approximations, the x and y components 
of the aircraft position are independently derived. Referring to 
FIG. 8, the value of the cross-track component x is determined by 
the position of the reference pattern image on the photosensitive 
array of the line-scan camera, the aircraft roll angle and the aircraft 
altitude as given by the following equation: 

x=z (tan .alpha.-.DELTA./f) where: 

z=true altitude above the runway; 

.alpha.=aircraft roll angle; 

.DELTA =displacement of the referenced pattem centerline image 
from the optic axis; 
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f=focal length of the camera lens. 

The distance .DELTA, is measured by multiplying the difference 
between the time intervals from the start of a line scan to the 
occurrence of the correlation peak, and that which corresponds to 
the reference pattern centerline image coinciding with the optic 
axis, by the line scan rate. 

The on-track component y of the aircraft position is measured by 
first determining from the laser aUimeter slant range and camera 
optics, the point P on the runway reference pattern where the size 
of its image on the camera's line array equals that of the 
correlator's reference pattern (see FIG. 3). At the instant the 
maximvim correlation occurs, the aircraft on-track position may be 
understood by reference to FIG. 9 and the following equation: 

y=y.sub.l +ztan .beta. 

where: 

y=aircraft distance to the shortest threshold point; 
y.sub.l ^distance from P to the threshold point; 
z=true altitude above runway; 
.beta =aircraft pitch angle. 

In the illustrated case, the pitch angle is negative as is typical with 
actual flight inspection patiis. Thus, the term z tan .beta, is also 
negative which results in y being smaller than y.sub.l. 

FIG. 10 illustrates the use of the optical pattern recognition system 
to update an inertial navigation system error state estimate based 
upon the precision position calculations of the optical pattern 
recognition system. 

FIG. 10 is a flow chart of the processing utilized to perform an 
update in a conventional Kahnan Filter. The optical pattern 
recognition system 102 (OPRS) provides the coordinates (x,y,z) 
where x is distance across runway centerline, y is distance along 
runway centerline and z is vertical. A standard coordinate 
conversion 104, in this case a simple rotation, converts these 
runway coordinates into precise east, north, and vertical 
coordinates (E, N, Z) .sub.REF of the aircraft at the instant of 
pattern passage. 

The INS 106 provides time tagged latitude (.phi.), longitude 
(.lambda.) and altitude (Z) which are interpolated to obtain the INS 
measured position at the time of pattern overflight. These geodetic 
coordinates are converted 108 into East, North and vertical 
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coordinates (E, N, Z) .sub.INS with respect to the plane tangent to 
the earth at runway pattern. Standard geodetic to planar coordinate 
conversions which account for the ellipticity of the earth are used 
for this step. 

At this point in the processing 1 10, we have two opinions of the 
position of the aircraft at the instant of pattern overflight, from the 
INS and the OPRS. The error in the INS measured position is 
typically on the order of hundreds of feet. The OPRS error is on 
the order of inches. The precision position obtained at each 
overflight of the runway pattern is used to make an update in the 
classic Kalman Fiher sense of the estimate of the errors in the 
Inertial Navigation System (INS). The simple algebraic difference 
between the two is an excellent estimate of the INS error for that 
instant of time. 

Therefore, this difference is used as the measurement vector in the 
standard Kalman Filter 1 12 on the error in the INS. In the Kalman 
Filter, the error state vector 1 14 is the estimate of the INS error 
states as to position, velocity, gyro drift rates, accelerometer 
biases, and mislevel. This error state vector is time varying and as 
a result, a time sequence of OPRS updates v^ll provide 
increasingly improved estimates of the states representing higher 
derivatives. The first pass provides excellent position, the second 
pass provides velocity and so on. 

These INS error state estimates are then used by the flight 
inspection computer to correct the INS measurements made from 
the Outer Marker dovra to the runway. These corrected INS 
measurements are used in tum as the reference for determination of 
the localizer and glidepath error in the Instrument Landing System 
being inspected. 

Also included in FIG. 6 is illumination system 76 which may be 
utilized to illuminate runway 10 and pattem 12 so that flight 
inspection can be accomplished during night hours. It should also 
be noted that in view of the fact that pattem 12 is not symmetrical, 
it will be necessary to provide an inbound/outbound signal to 
correlator 66 so that the reference pattem template may be reversed 
to detect pattem 12. 

With reference now to FIG. 7, there is depicted a more detailed 
block diagram of the video position reference system of FIG. 5. As 
can be seen, digital line scaiming camera 24 produces, in the 
depicted embodiment, two output video channels 76 and 78. In a 
preferred embodiment of the present invention, digital scanning 
camera 24 is implemented utilizing an industrial line scan camera 
such as the Model CCD 1300R manufactured by Fairchild Camera 
and Instrument Corporation of Palo Alto, Calif. This particular 
model utilizes CCD imaging techniques and puts out two time 
division multiplex video outputs. 
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In order to extract the essential features of geometric pattern 12 
over a wide range of absolute image brightness, analog video 
outputs 76 and 78 are coupled into a pair of wide-band logarithmic 
response amplifiers 84 and 86. In a preferred embodiment of the 
present invention an independent auto-biasing circuit is utilized in 
conjxmction with each logarithmic response amplifier to establish a 
precise black level operating point. This is achieved by activating 
the sample state of sample-and-hold circuits 80 and 82 at a pixel 
count which coincides with the black reference video output level. 

The outputs of logarithmic response amplifiers 84 and 86 are then 
coupled to linear difference amplifier 88. Those ordinarily skilled 
in the art will appreciate that the well known mathematical 
properties of logarithms will result in the output of difference 
amplifier 88 being a constant value for sets of video input signals 
having a constant ratio, regardless of absolute value. Since the two 
video channel odd/even format produces light intensity levels of 
adjacent pixels simultaneously, the difference amplifier output 
indicates the ratio of light intensity change between adjacent 
pixels. However, a polarity reversal is introduced each pixel period 
because of the altemating sign of the most recent data. 

Although a sharp black-to-white or white-to-black transition is the 
video feature being sought, allowance must be made for the 
possibility of the edge of geometric pattern 12 falling in the center 
of a pixel sensing element producing an intermediate average 
intensity. For this reason, the ratio threshold value must be set at 
somewhat less than a ratio of two-to-one to ensure adequate 
probability of detection. Voltages corresponding to the threshold 
ratio and its reciprocal are supplied to a pair of analog comparator 
amplifiers 90 and 92 which are coupled to the output of high pass 
filter 94. The outputs of analog comparator amplifiers 90 and 92 
are then decommutated on an odd/even basis utilizing switch 96. 
The resultant ouptut enables the setting and clearing of bistable 
multivibrator 98 which produces a one-bit quantized video output 
which is required by digital correlator 66. In the depicted mode of 
the present invention, the state of that bit will be zero for a black 
detection and one for a white detection. 

As explained above, a digital reference pattern is stored within 
reference pattem register 100 and then correlated with the digitized 
output of a line scanning camera to determine the detection of 
geometric pattem 12. Digital peak detector 68 is implemented 
utilizing a seven-bit latch and is updated each time the output of 
correlator 66 is greater than a previous output. Oscillator 74 is 
utilized to generate timing pulses, and cross-track register 70 and 
scan number register 72 are implemented utilizing an eleven-bit 
latch and a four-bit latch respectively. Eleven-bit counter 71 is 
utilized in conjunction with cross-track register 70 and four-bit 
counter 73 is utilized in conjunction with scan number register 72. 
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Those ordinarily skilled in the art will appreciate upon reference to 
the foregoing specification that when utilizing the video position 
reference system of the present invention of an inspection aircraft 
may fly over an airport runway and utilize visual information to 
obtain a substantially exact reference position of the aircraft with 
respect to the runway which may then be utilized to upgrade the 
onboard inertial navigation system to an accuracy sufficient to 
permit Class III instrument landing systems to be calibrated. 

Although the invention has been described with reference to a 
specific embodiment, this description is not meant to be construed 
in a limitmg sense. Various modifications of the disclosed 
embodiment as well as alternate embodiments of the invention will 
become apparent to persons skilled in the art upon reference to the 
description of the invention. It is therefore contemplated that the 
appended claims will cover any such modifications or 
embodiments that fall within the true scope of the invention. 
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Observational Evidence for an Active Surface Reservoir 
of Solid Carbon Dioxide on Mars 

Michael C. Malin,* Michael A. Capllnger, Scott D. Davis 

IHigh-resolution images of the south polar residual cap of Mars acquired in 1999 and 
2001 show changes in the configuration of pits, intervening ridges, and isolated 
mounds. Escarpments have retreated 1 to 3 meters in 1 martian year, changes that 
are an order of magnitude larger than can be explained by the sublimation of water 
ice, but close to what is expected for sublimation of carbon dioxide ice. These 
observations support a 35-year-old conjecture that Mars has a large surface 
reservoir of solid carbon dioxide. The erosion implies that this reservoir is not in 
equilibrium with the present environment and that global climate change is occuning 
on Mars. 

Malin Space Science Systems. Post Office Box 910148, San Diego, OA 92191- 
0148, USA. 
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In their seminal 1966 paper on the behavior of CO2 and other volatiles on Mars, 

Leighton and Murray (1) deduced from a numerical thermal model that the polar 
caps are composed of frozen CO2. They drew three other important conclusions 
from this analysis: (i) The year-round presence of solid CO2 regulates the pressure 

of the atmosphere, (ti) The atmospheric pressure changes cyclically on a semiannual 
basis because of fomiation and retreat of the seasonal frost caps. (\\\) The total 
amount of CO2 on Mars could greatly exceed that presently in the atmosphere and 

form a large surface resen/oir. in the past 35 years, many of Leighton and Murray's 
predictions have been verified and demonstrated by measurements taken on Mars 
(2). However, the existence of a surface reservoir of CO2 large enough to have 

important implications for long-term climate and climate stability has yet to be 
demonstrated (3, 4). Here we report observations, made with the Mars Global 


http .//www. sciencemag. org/cgi/content/full/sci ;294/5 549/2 1 46 


5/24/07 


Surveyor (MGS)/Mars Orbiter Camera (MOC), of interannual morphologic changes 
in the south polar residual cap that support the argument for a dimatologically 
important CO2 reservoir. 

The south polar residual cap exhibits morphologies unique to that location, 
characterized by irregular to circular pits, remnant mesas, and other landforms, 
interpreted to have formed by collapse and erosion (5). These features range from a 
few tens of meters to a few hundreds of meters in scale but display only a few 
meters of relief (5). 

The erosion appears confined to, and highlights the layered nature of, the uppermost 
materials that make up the cap. The layers are relatively uniform in thickness but 
display a range in surfece expression that appears related to their ability to erode 
(6). As few as one to as many as six layers crop out at various locations within the 
cap; not all layers are seen at all locales. On the basis of shadow measurements 
made during this investigation, each layer is about 3 m thick. 

MOC was first used to observe the south polar region at high resolution in July 
1999 (7). To search for changes, specific locations within the residual cap were 
reimaged at the same season and time of day beginning in late July 2001 [starting 
around heliocentric longitude (L^) 223°±4° and continuing since then]. In each case, 

a 2.2 to 3.7 nfi/pixel image covering a portion of an earlier image was acquired (8). 
Upon receipt, both the original and new images were processed for analysis (9). 

Figure 1 shows good examples of pits, intervening ridges, and isolated mounds that 
have changed. Diametric measurements, made on 100 features on each of four 
image pairs, indicate that these pits enlarged or their intervening ridges shrank by 
about 6 ± 2 m and features smaller than -6 m across disappeared. In other words, in 
these examples, nonhorizontal surfaces appear to have retreated by ~3 m. Other 
areas within the residual cap display different morphologies, and some of these 
show no detectable change. Preliminary statistics suggest that between 25 and 50% 
of the escarpments have retreated between 1 and 3 m. 


Fig, 1. Pits, intervening ridges, and mounds in south polar 
residual ice cap. Extracted from orthographically projected, 
1 m/pixel versions of M08-04516 (left; 10/19/1999, ~ 227^) 

and E07-01 565 (right; 8/23/2001 , L3 21 9"). Pit walls and other 

escarpments have retreated 3 m in one Mars year. 
Illumination from upper left. (View Larger Version of this Image 
(84K GIF file)] 
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As part of the input to an estimate of the total amount of scarp retreat that occurred 
during the past martian year, the total perimeter in a variety of settings was 
detemiined by image processing techniques (10) and divided by the area viewed. 
Measurements were made on images representative of the range in areal density of 
various forms of escarpments within the residual south polar cap. The average scarp 
perimeter per area was 2.4 x 10"^ m/m~^. 


http://www.sciencemag.org/cgi/content/full/sci;294/5549/2146 


5/24/07 


We determined that the amount of retreat was consistent with the subliming of a 
volatile ice by calculating the annual solar insolation using software based on 
spacecraft navigation routines, and we compared our results with the results of 
previous models of polar cap behavior (1, 11-13 . 14). We used the order of 
magnitude greater volatility of CO2 relative to water ice to differentiate between these 

two candidates. Table 1 summarizes the results of our calculations of the insolation 
and the equivalent areal mass density for H2O and CO2 ice for horizontal surfaces at 

various latitudes from the equator to near the south pole. 


Table 1. Calculated solar insolation and maximum equivalent ice mass erosion 
for HjO ice (column 6) and CO2 ice (column 7) for one martian year on horizontal 

surfaces at different latitudes (column 1). Columns 2, 3, 4, and 5 give computed 
solar insolation in energy areal density per Mars year, or its equivalent power 
density, in several different representative physical units. All four columns present 
the same information-the solar energy input into the surface, independent of 
surface material-but in different units. Column 6 was computed using 652 cal/g 
for the latent heat of H2O at 145 K. Column 7 was computed using 143 cal/g for 

the latent heat of CO2 at 145 K. Values from columns 6 and 7 are used as input to 

the calculations summarized in Table 2. 


Parameter 


Maximum 








CO2 

Latitude 

(joules m ^ 

(Wm 

(cal m ^ 

(cal cm ^ 

(g cm~2 

(g cm~2 


year~^) 


year~^) 

year~^) 

year~^) 

year"^) 

-85° 

4.79 X 10^ 

80.7 

1.14 X 10^ 
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4.89 X 10^ 
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1.07 X 10^° 

179.6 

2.54 X 109 

2.54 X io5 

389 

1775 


Assistant Chair, 
Department of 
Neurobiology and 
Physiology, Northwestern 

University 

Northwestern University 
Chicago, IL 

Chair of Developmental 
Biology and Anatomy 
University of Edinburgh 
Edinburgh, United Kingdom 

Great Lakes Invasive 
Species Research Scientist 
Central Michigan University 
Mount Pleasant, MI 

Cancer Translational 

Scientists 

AstraZeneca 

Alderley Park, Cheshire, 

South Manchester, United 

Kingdom 

Research and Development 
Scientist 

Saladax Biomedical Inc 
Bethlehem, PA 

DIRECTOR FOR RESEARCH 
University of Texas Health 
Science Center, San 
Antonio 

San Antonio, TX 
CHAIR 

Duke University 
Durham, NC 

More jobs 


A variety of factors affect the translation of the maximum equivalent areal mass loss 
(Table 1 ) to a comparable effective scarp retreat (Table 2 ). and the result is sensitive 
to the values assumed for each of these factors. We explored some of these factors 
and a range of possible values for them (albedos between 0.5 and 0.75 and 
densities between 1 and 3 g/cm^). Atmospheric contributions (which can act to 
enhance or diminish sublimation, depending on seasonal factors), conduction into 
the surface, and other physical processes that would reduce the amount of energy 
available for sublimation were represented by an efficiency factor of arbitrary value. 
We also assumed that slopes of between 30** and 50"* (typical of the slopes 
measured by shadow techniques) facing the equator at SS^'S latitude receive roughly 
the same insolation (with the same efficiency) as horizontal surfaces at 35'' to SS^'S 
(this simplification encompasses variations resulting fi^om slope azimuth and 
season). 


Table 2. Calculated amount of ice mass erosion per Mars year for and GO^ 
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ice as a function of latitude (column 1) and energy transfer. Columns 2 and 

3 show the amount of mass loss per unit area for an albedo of 0.65. Columns 

4 through 9 show the equivalent amount of surface retreat for material with a 
density of 1 .5 g/cm and three different efficiencies of energy transfer 

(100, 50, and 80%). Calculations are for horizontal surfaces; escarpment retreat 
can be estimated by subtracting the slope declivity for the "actual" latitude to 
derive an "effective" latitude and reading across the row of that latitude. 

Parameter 

Efficiency 

= 100% 

Efficiency 

= 50% 

Efficiency = 80% 


CO2 

HjO 

CO2 

H2O 

CO2 

HjO CO2 
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(g cm~2 

(cm 
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(cm (cm 

year"^) 

year"^) 

year~^) 

year"^) 

year ^) year ^) 

year~^) year"^) 

-85° 61 

279 

41 

186 

20 

93 

33 149 

-80° 62 

285 

42 

190 

21 

95 

33 152 

-75° 65 

295 

43 

196 

22 

98 

34 157 

-70° 68 

309 

45 

206 

23 

103 

36 165 

-45° 101 

462 

68 

308 

34 

154 

54 246 

0° 136 

621 

91 

414 

45 

207 

73 331 


We first compared our calculations for horizontal surfaces with the areal mass 
density of CO2 frost computed using thermal models of the seasonal polar cap near 

the south pole. From these models, we calculated that between 100 and 150 g/cm^ 
is deposited at 80''S latitude each winter and is removed by sublimation each spring 
and summer (1, 14, H, 12). When we use parameters similar to those used in these 
models (albedo = 0.65, density = 1.0 g/cm^), efficiencies between 35 and 50% 
reproduce their results; an efficiency of less than 35% permits net accumulation at 
SO^'S, whereas a greater efficiency permits net erosion or the heating of the 
nonvolatile surface. 

For the same assumptions, our calculations suggest that efficiencies between 
40 and 80% will lead to 1 to 3 m of retreat of 30"* to 50'' escarpments in CO2. More 

importantly, they show that scarps in water ice could, at most, retreat only a few 
decimeters. Only extreme combinations of our three parameters (albedo < 0.5, 
density < 0.5 g/cm^, and "efficiency" > 80%), which are generally unrealistic, could 
reproduce in water ice the magnitude of scarp retreat actually observed. We 
conclude that the erosion is occurring in one or more layers of moderately dense 
CO2 ice. 

We estimated the amount of material lost from the south polar cap over the past 
year by mapping the extent of the pitted, layered material (5) on a south polar 
stereographic map projection of a Mars digital image mosaic, based on close 
inspection of 144 pre- and immediately postsunrise, 12 to 15 m/pixel survey images 
acquired during July and August 1999. The map was then transferred to a 
comparably projected late summer view of the south polar region acquired by the 
MOC wide angle camera. The extremely strong correlation of the mapped contacts 
with the albedo of the late summer residual cap in the wide angle image permitted 
the use of image processing techniques to determine the area of the cap and its 
pitted surface (-87,000 km^). The calculated amount of scarp perimeter is -2x10^ 
m, and the volume of material lost between 1999 and 2001 is -2 to 10 x 10^ m^ (2 to 
10 X 10^5 cm^) for the ranges in retreat distance (1 to 3 m) and fraction of 
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escarpments experiencing retreat (25 to 50%). For densities of COg ice of 1 to 

2 g/cm^, this is 2 to 20 x 10^^ g, or about 0.008% to 0.08% of the mass of the 
atmosphere. 

A 3 m thick layer of pure COj ice with a density between 1 and 2 g/cm^ covering the 
entire residual cap would have a mass of 2 to 4 x 10^^ g (1 to 2% of the present 
mass of the martian atmosphere). At a linear recession rate, retreat of the present 
escarpments would consume a layer in 30 to 150 martian years. However, 
measurements of the enlargement of different sized pits show that the retreat of 
escarpments is independent of the size of the pit, so the rate of pit growth and hence 
release of additional CO2 into the atmosphere is not linear but will increase with time. 

For example, the perimeter of a 10*m-diameter pit would double in 2 martian years 
at a 3 m/year recession rate and at that point would experience twice as much 
erosion and contribute twice as much gas to the atmosphere as It did 2 years earlier. 
The size-frequency and areal density relationships for depressions within the polar 
cap have not been fully quantified, but preliminary values suggest that the erosion of 
a layer may take only one-third as long as the linear estimate. 

Thus, without accounting for deposition elsewhere, this erosion rate implies a 
secular and increasing contribution of C02into the atmosphere, augmenting its 
mass by as much as 1% per martian decade. Were sufficient layers of CO2 

available, it would only take a few martian centuries to double the present 
atmosphere. 

We do not know where the COj is going (15). As previously noted (3). 
thermodynamics favor precipitation of CO2 in the north polar region, but CO2 frost is 
not seen to persist on the northern cap through the summer. The recently released 
CO2 probably participates in the formation and retreat of the seasonal frost caps. 

Where it goes from there is uncertain; as the atmospheric pressure rises, regolith 
adsorption may increase (16, 17) or COj may be precipitated at locations 

somewhere in the polar regions not presently detected. 

The large, cyclic, and potentially chaotic obliquity and orbital eccentricity variations 
that Mars is believed to experience as a result of gravitational interactions between 
the planets may play a role in how, where, and when the surface reservoir forms, but 
calculations (18-22) show that values for these attributes have not changed 
substantially over the past 10,000 to 20,000 years. Given the high rates of 
destruction of the CO2 layer(s), it is not clear that any of them could be surviving 
from the last period of dramatically low obliquity (when the polar regions would have 
been substantially colder) or dramatically high obliquity (when the higher polar 
temperatures might have been offset by higher atmospheric pressures to permit 
polar deposition) or from the last period when the polar temperatures were changing 
rapidly. 

Repeated MOC observations indicate that a surface reservoir of solid CO2 persists 

year-round on Mars and that this deposit is not in equilibrium with the present polar 
environment (it is being eroded). The rate and amount of erosion are phenomenal- 
the amount of scarp retreat, although small on a planetary scale, is substantially 
greater than any change previously identified on Mars at meter scale (7). It is 
unlikely that this rate can be sustained indefinitely, because the entire relief beneath 
the visible portion of the south polar residual cap, were it CO2, could be sublimed to 
vapor in a few thousand Mars years. These and other observations (7, 23) suggest 
that the present martian environment is neither stable nor typical of the past. 
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